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the mean king problem 



o 
o 

(N 

> 
O 



Juan Pablo Paz, 1,2 Augusto Jose Roncaglia, 1, 2 and Marcos Saraceno 3 ' 4 

1 Departamento de Fisica, FCEyN, UBA, Pabellon 1, 
Ciudad Universitaria, 1428 Buenos Aires, Argentina 
^Theoretical Division, LANL, MSB213, Los Alamos, NM 87545, USA 
Unidad de Actividad Fisica, Tandar, CNEA, Buenos Aires, Argentina 
^Escuela de Ciencia y Tecnologia, Univ. Nac. Gral. San Martin, Buenos Aires, Argentina 

(Dated: February 1, 2008) 

We analyze and further develop a new method to represent the quantum state of a system of n 
qubits in a phase space grid of N x N points (where N = 2 n ). Th e method, which was recently 
proposed by Wootters and co-workers (Gibbons et ah, quant-ph/0401155 ), is based on the use of 
the elements of the finite field GF(2 n ) to label the phase space axes. We present a self-contained 
overview of the method, we give new insights on some of its features and we apply it to investigate 
problems which are of interest for quantum information theory: We analyze the phase space repre- 
sentation of stabilizer states and quantum error correction codes and present a phase space solution 
to the so-called "mean king problem" . 

PACS numbers: 3.65.Ca, 3.65.Ta, 3.65. Wj, 02.10.De 



(N 
> 



o ; 
^ ■ 

Oh! 
> ■ 

C ■ 

=5 

cr 



X 



I. INTRODUCTION 

Quantum mechanics can be formulated in phase space, 
the natural arena of classical physics. When doing this, 
quantum states are represented by phase space distribu- 
tions with peculiar properties that distinguish them from 
their classical counterparts. The most common one is the 
Wigner function pj which is a real-valued function that, 
unlike a genuine probability density, takes negative val- 
ues for generic quantum states. However, Wigner func- 
tions share an important property with a classical prob- 
ability density: when integrated along any phase space 
strip (the region between any two parallel lines) is a pos- 
itive number bounded between zero and one. The use of 
phase space methods has proved useful in the study of 
the quantum-classical transition and in the analysis of 
semi-classical properties of quantum systems. 

In the context of the recent surge of interest on quan- 
tum information the use of phase space representation 
for quantum states and algorithms has been proposed 
0. However, as opposed to what happens for contin- 
uous systems, the definition of Wigner functions for the 
discrete case suffers from some ambiguities. In fact, there 
are various approaches available to generalize the Wigner 
function for quantum systems with a finite-dimensional 
space of states. Discrete versions of Wigner functions 
were first introduced in the context of studies of semi- 
classical properties of classically chaotic systems by Han- 
nay and Berry [3|]. Some time later, Feynman |4| indepen- 
dently defined a phase space distribution for a spin-1/2 
particle in a 2 x 2 grid (see also ||). Wootters fur- 
ther developed this idea and defined the discrete version 
of the Wigner function for JV-dimensional quantum sys- 
tems when N is a prime number (see also 7]). In such 
case the phase space is an N x N grid. Wootters method 
can also be applied when N is a composite number but in 
such case the Wigner function must be defined in a phase 



space grid which is the Cartesian product of the phase 
spaces corresponding to the prime factors of N. For the 
prime dimensional case, the phase space points (q, p) have 
position and momentum coordinates which are integers 
between and N — 1. The fact that this set has the 
structure of a finite field is essential in implementing the 
method. A different approach was followed by Cohendet 
et al. 8] who were able to define a Wigner representa- 
tion for odd values of N (not necessarily prime) using an 
N x N grid. Similarly, Leonhard [j| extended the idea 
for even values of N but in this case it was necessary to 
use a grid with 2N x 2A*" points. Some of us proposed 
a similar approach 2] which can be used to define a dis- 
crete Wigner function for arbitrary values of N using a 
grid with 2N x 2N points. In all these cases the phase 
space coordinates (q,p) are assumed to be integers be- 
tween and N — 1 (or 2N — 1). Arithmetic modulo N 
is used when operating with such coordinates. The fact 
that for non-prime values of N Zm is not a finite field 
is the cause of several peculiar and unpleasant properties 
of the phase space (for example, when defining lines in 
phase space as solutions of linear equations one ends up 
with non-parallel lines that intersect in more than one 
point). More recently Wootters [Tjj proposed a differ- 
ent approach to define the discrete Wigner function that 
enables one to work on a phase space grid with N x N 
points when N is the power of a prime number. The 
method, which was discussed in a comprehensive paper 
by Gibbons, Hoffman and Wootters [llj, heavily relies 
on the use of the elements of the finite field GF(N) to 
label both phase space coordinates (q,p). For various 
reasons this approach seems to be well suited to study 
some quantum information problems. Thus, quantum 
computers made out of n qubits have a Hilbert space 
whose dimension N = 2™ is a power of a prime num- 
ber. Moreover, operators made out of tensor products 
of Pauli operators acting on each qubit play an impor- 
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tant role for quantum computers and have a central role 
in Wootters phase space method (they represent phase 
space translation operators as discussed below). 

The use of the discrete Wigner functions in quantum 
information has been recently been proposed. Applica- 
tions using Wigner functions in a 2N x 2N grid include: 
the study of general properties of quantum states and the 
phase space representation of Grover's search algorithm 
|2|, the phase space representation of quantum telepor- 
tation of arbitrary dimensional systems [12j, the study 
of quantum walk algorithms (where phase space meth- 
ods seem to be particularly well suited) |13|. the anal- 
ysis of decoherence models with a natural phase space 
representation |l4j . applications to the use of quantum 
computers to study properties of classically chaotic maps 
[TH and the development of efficient techniques for quan- 
tum state tomography in phase space 0, . On the 
other hand, some attention was devoted to the use of 
the Wigner function o rigi nally defined by Wootters [a] to 
analyze teleportation jl8| . Following Wootters's more re- 
cent ideas [lij a study of the nature of the set of quantum 
states that have positive Wigner functions was presented 
by Galvao who also conjectured the existence of a 
connection between these states and the ones that can 
be classically simulatable. However, it is still unclear if 
the use of phase space methods will help in expanding 
our understanding of some of the problems of quantum 
information. It is indeed possible that the various def- 
initions of Wigner function may end up proving to be 
useful to analyze different problems (indeed, much of the 
use of 2N x 2N phase space representations was done in 
the spirit of previous studies of semi-classical properties 
of quantum maps [3j). 

In this work we will present a new step towards ex- 
ploring the use of phase space techniques for quantum 
information using and expandin g th e recent work done 
by Wootters and co-workers 0, lid ITl| . As mentioned 
above, this method relies heavily on the fact that phase 
space coordinates (position and momentum) are chosen 
to be elements of a finite Galois field. In this way it is 
simple to construct a phase space with geometric proper- 
ties that are similar to the usual phase space (i.e., where 
two lines that are not parallel have only one intersection, 
etc). Once the phase space arena is built in this way it is 
possible to define ways to map sets of parallel lines (i.e., 
striations) to bases of the Hilbert space, and to show 
that the bases associated with different striations turn 
out to be mutually unbiased. In our paper we will dis- 
cuss some interesting properties of this phase space for 
systems made out of qubits. We will analyse a few prob- 
lems where these phase space tools could be naturally 
applicable. 

The paper is organized as follows. In Section II we re- 
view, in a self-contained way, the main properties of the 
discrete phase space. Here, we also show how to reduce 
the arbitrariness in the association between lines and 
states (fixing the quantum net) by imposing additional 
symmetries to the phase space structure. In Section III 



we show how to define the discrete Wigner function and 
discuss some of its most important properties. In Section 
IV we analyze quantum information problems using this 
discrete Wigner function: First we analyze properties of 
stabilizer states |20j , which are eigenstates of translation 
operators. Finally, we show a phase space solution to 
the mean king problem [2lLl22j|. In Section V we analyze 
some general properties of the phase space representation 
of stabilizer states, we discuss possible future directions 
and summarize our results. 



II. DISCRETE PHASE SPACE FOR A SYSTEM 
OF n QUBITS 

A. Phase space coordinates 

In this section we will present a self contained review 
of a recent method to define a phase space for a system of 
n qubits. This discrete phase space is the natural arena 
in which a discrete Wigner function can be introduced. 
Our treatment is based on the ideas discussed first by 
Wootters ly] and developed by Gibbons, Hoffman and 
Wootters [Tot ITU ]. We will concentrate on systems of 
qubits but the method can be extended to systems with 
Hilbert space with a dimension that is a power of a prime 
number. 

For a system with n qubits we introduce a phase space 
grid with N x N points (N — 2 n denotes the dimension 
of the Hilbert space). We label each phase space point 
with its "position" and "momentum" coordinates (q,p), 
where both q and p are elements of the field GF(2 n ). 
The elements of this field can be thought of as being 
the N = 2 n n-tuples with binary (i.e., or 1) entries. 
The sum in GF(2 n ) is defined as the bitwise (mod 2) 
addition of the respective n-tuples. Elements of GF(2 n ) 
can also be thought of as polynomials of degree n — 1 
with binary coefficients. The product in GF(2 n ) is de- 
fined as the product of the corresponding polynomials 
modulo a primitive polynomial (this is a polynomial of 
degree n which is irreducible, i.e. it cannot be factored 
as the product of two polynomials of lower degree). The 
product in GF(2 n ) is such that the set formed by all 
the elements of GF{2 n ) excluding the zero is a cyclic 
group of order 2™ — 1 with the multiplication in GF(2 n ) 
as the group operation. Thus, GF(2 n ) contains the zero 
element and the 2™ — 1 different powers of a primitive cl- 
ement lj: {0, 1, uj, ui 2 , . . . , a; 2 "" 2 } (note that uj 2 "^ 1 = 1). 
This establishes a natural order for the field elements, 
and we use this order to label the axes. Using the primi- 
tive polynomial it is clear that any power ui k with k > n 
can be expressed as a linear combination of powers of 
degree lower than n. This defines the addition rule in 
the field. For example, for n — 3 (three qubits) the 
polynomial tt(x) — x 3 + x 2 + 1 can be chosen as primi- 
tive. Then, the field is formed by the following elements 
GF(2 3 ) = {0,l,w,w 2 ,l+w 2 ,l+w + w 2 ,l+w,w + w 2 } 
(the last four elements in the list are the powers cu k with 
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k = 3,4, 5, 6). An ordered set of n elements of the field, 
E = {eo, . . . , e„_i}, is a basis if every element of GF(2 n ) 
can be written as the linear combination 

n— l 

x = Xjej. (1) 

i=0 

In this equation all the coefhcients elements of 

GF(2) (i.e., they are either or 1). The dual basis of 
E can be introduced as follows: The trace of any field 
element x is defined as 

tr(ar) = x + x 2 + x 2 ' + . . . + x 2 "~\ (2) 

This operation is linear, which follows from the fact that 
for every pair of elements of GF(2 n ) the square of their 
sum is equal to the sum of their squares: i.e., [x\ +X2) 2 — 
x\ + x\. Moreover, it can be shown that the trace (0) 
takes values in GF{2). Using the above trace one can 
define a dual for every field basis: Thus, it turns out that 
for every basis E there is a unique basis E satisfying that 
tr(eiej) = 6ij. The basis E is the dual of E. 

There is a very useful representation for GF(2 n ) in 
terms oftixn matrices with binary coefficients. It is 
worth pointing out here the main properties of this rep- 
resentation, which are discussed in more detail in the Ap- 
pendix. Given a primitive polynomial tt(x) for the field 
GF(2 n ) one can choose its companion matrix M as the 
primitive element for the matrix representation (such ma- 
trix is defined in the Appendix and is such that tt(M) = 
and M 2 _1 = 1). The first n powers of the matrix M can 
be chosen as the canonical basis for the field and any field 
element is then written as x = X)"=o XiM 1 with the n 
binary coefficients Xi being the coordinates representing 
x as an n-tuple. The virtue of the matrix representation 
of GF(2 n ) is that both the sum and the product in the 
field correspond directly to ordinary operations between 
matrices. Also, the trace of any element, defined in J2J is 
nothing but the trace of the corresponding matrix. These 
are the basic mathematical elements about GF(2 n ) that 
are necessary to build the phase space arena to represent 
the quantum states of n qubits. 



0)" 


v 


6 


5 


4 


3 


2 


1 





CO ' 


V 


5 


4 


3 


2 


1 





6 


CO 4 


V 


4 


3 


2 


1 





6 


5 


CO 3 


V 


3 


2 


1 





6 


5 


4 


CO : 


V 


2 


1 





6 


5 


4 


3 


CO 


V 


1 





6 


5 


4 


3 


2 


1 


V 





6 


5 


4 


3 


2 


1 





■ 




h 


h 


h 


h 


h 


h 







1 


CO 


CO 2 


CO 3 


CO 4 


CO 5 


CO" 



FIG. 1: Lines passing through the origin (rays) of phase space 
for 3 qubits. All rays apart from the vertical and horizontal 
ones (labeled with h and v in the plot) satisfy the equation 
p = u 3 q (the power j defines the slope of each curve and its 
value is displayed inside each cell) . The total number of lines 
in this case is N + 1 = 9. The position and momentum axes 
are labeled by elements of the field GF(2 3 ). 

Also, because N is a power of a prime number, two lines 
that are not parallel intersect only at one point. 

A phase space translation by an amount ctQ — {qo,Po) 
will transform any phase space point a = (q,p) into 
TagCt — a + ao- Lines in phase space are invariant un- 
der some translations. Consider the line containing the 
origin and the point (q,p). This line is formed by the 
points (sq,sp), where s is a parameter that ranges over 
all the elements of GF(2 n ). The striation containing this 
line remains invariant under any translation of the form 
T (sq,sp)- We will take as a reference line for each stria- 
tion the line passing through the origin, specified by the 
equation aq + bp = 0, which we call a ray Notice 
that if the axes are labeled with the powers of a primi- 
tive element these rays appear on the main diagonals and 
look "parallel" in the ordinary sense, wrapping periodi- 
cally on a torus of 2™ — 1 periodicity. The N + 1 rays for 
three qubits are displayed in FigQ 



C. Translation operators 



B. Lines and translations in phase space 

In the N x N phase space grid we define a line as the 
set of phase space points satisfying a linear equation aq + 
bp = c (where all elements and operations appearing in 
this equation belong to GF(2 n )). Two lines are parallel 
if they do not intersect, or equivalently if they satisfy the 
equations aq + bp = c and aq + bp — d with c^c'. Lines 
in this phase space grid have rather natural properties: 
Given a line it is always possible to find N — 1 other lines 
that are parallel to it. A set of N parallel lines form a 
striation of the phase space (i.e., a striation is a foliation 
of the grid with parallel lines). It is easy to show that 
there are N + 1 different striations of the phase space. 



The association of a translation operator T a with a 
phase space point a = (q,p) should be done in such a 
way as to preserve the additive structure of the field. 
Moreover, following we require that for a sys- 

tem of n qubits translations should act independently on 
each qubit, thus preserving the tensor product structure 
of the Hilbert space (this is a nontrivial assumption which 
is not satisfied in other constructions |2j). At the level 
of individual qubits position and momentum translations 
are identified with Pauli operators Xi and Zi. Operators 
that satisfy these requirements are 

n-l 

^(q, P) = II XitoZi**^**-** = A q Z p e^ q ' p , (3) 

i=0 
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where we denote (q, p) the binary strings 
{qo ■ ■ -q n -i,Po ■ ■ -Pn-i)- The phase q p = J^liPi 
is added to make the operators both unitary and 
hermitian. It is simple to show that the factor 
exp(z7rq • p/2) can take values ±1, +i and that 
T(q,0) = X q , T(0,p) = ZP, T(q,q) = F q . The set of 
iV 2 operators T(a, b) is an orthogonal basis of the space 
of operators since Tr(T(a, b)T(c, d)) = NS^^S^.d)- 

Associating a translation T(q, p) with a phase space 
point a = (q,p) requires a mapping between the field 
elements q and p and the binary strings q and p. The 
most natural such mapping is to use the binary string 
q (p) formed by the coefficients of the expansion of the 
field elements q (p) in a given field basis. The field basis is 
arbitrary and could in principle be different for position 
and momentum. However, we will show below that the 
consistency of the phase space structure implies that once 
we choose a basis to expand the position q, there are 
strong constraints on the basis we could use to expand 
the momentum p. Before doing this, it is worth noticing 
two simple properties of the translation operators defined 
above: The composition law for the translations is 

T(a, b)T(q, p) = ±T(a + q, b + p) exp(i~(b ■ q- a- p)), 
which leads to the commutator 

[T(a, b) , T(q, p)] = ±2i sin - (b • q - a • p)T(a + q, b + p) . 

(5) 

Thus two translations commute iff a • p — b • q = 
(mod 2). 

Using the results in the Appendix we can then con- 
struct commuting sets of translations as follows: consider 
the set of translations 

T(aM j , bM j ) j = 0, 2 n - 2, (6) 

where M is the companion matrix of the primitive poly- 
nomial used to construct the field and M is its transpose. 
Using Eq. (J5J it is easy to show that these operators 
commute. These operators, together with T(0, 0) form 
a complete set of N commuting translations whose com- 
mon set of eigenvectors constitute a basis for the Hilbert 
space. Using this procedure we can partition the set of 
all N 2 translation operators into N + 1 sub-sets of com- 
muting translations. Each of these sub sets contains the 
identity and other N — 1 operators which are of the form 
©. Thus, each of these N + 1 sub-sets can be asso- 
ciated with different strings (a, b). The following is a 
simple way to choose the N + 1 strings to define the cor- 
responding sets: We can first define two sets associated 
with the strings (1,0) and (0,1). Then, the remaining 
N — 1 sets can be associated with strings of the form 
(l,b) for all b ^ 0. Here, and below, we use the nota- 
tion 1 = (10 ... 0) and = (00 ... 0). We remark that 
this is a very simple way to find the N+l commuting sets 
of translations defining the corresponding MUBs. The 
above choice of binary strings a and b is arbitrary but it 



is easy to see that different choices of such strings would 
lead to the same N + l sets. There are two extra free- 
doms that one may exploit. First, one can change the 
companion matrix M by choosing a different primitive 
polynomial. Second, one can define the translation oper- 
ators in a different way than the one given in eq. © . For 
example, for each qubit we could interchange the three 
Pauli operators. In this way one obtains a different par- 
tition of the Pauli group into N+l commuting sets. 

D. Association between lines and states 

With the above tools we are ready to address the first 
fundamental point made by Wootters in defining the 
phase space. The key to this construction is to estab- 
lish a one to one correspondence between every phase 
space line A and a state in Hilbert space. We will de- 
scribe how to associate a pure state with the line A (or, 
analogously, define a mapping between every line A and a 
rank one projection operator P(\)). Following Wootters 
this mapping is defined by imposing a natural geometric 
constraint: we require that the mapping P(X) should be 
such that the state associated with the translated line 
T a A should be obtained from P(X) by applying a trans- 
lation operator T a , i.e. 

P(r a X) = T a P{\)Tl (7) 

This condition enforces the validity of two important 
results. The first result following from 10 is that lines 
belonging to the same striation are associated with or- 
thogonal states (i.e., a striation is associated with an or- 
thonormal basis of the Hilbert space). The reason why 
this is the case is the following: Consider the striation 
containing the ray formed by the points {sa, sb) where a 
and b are fixed and s ranges over all elements of GF(2 n ). 
All lines in this striation are invariant under translations 
r a with a = (ta, tb),Vt € GF(2 n ). Therefore, as r a X = A, 
equation (J7J implies that all the operators T a must com- 
mute with the projectors P(X). Thus, this implies that 
the translation operators associated with every point in 
the ray (sa, sb) must form a commuting set and, more- 
over, that the states associated with the striation con- 
taining the ray must be the common eigenstates of the 
translations. 

The second result following from Q is that the bases 
associated with the iV+1 different striations are mutually 
unbiased. In fact, one of the virtues of the phase space 
structure based on finite fields is that it establishes a clear 
connection between phase space striations and mutually 
unbiased bases (MUB). This concept is interesting on its 
own and has been studied in detail [2^, 0, H^, |26j . Two 
bases sets {\fa), j = 1, • • • , N} and {\i(>k), k = 1,.. ., N} 
are MUB if and only if | {fa \tp k ) | 2 = 1/N for all values of 
j and k. It has been shown that N + l MUB exist when 
N is the power of a prime number |24|. The interest 
in MUB is strongly connected with the problem of state 
tomography. Thus, one can show that the most efficient 
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way of completely determining a quantum state is by 
making von Neumann measurements on the N + 1 MUB 
[24l | . Notably, Wootters phase space construction is itself 
a relatively simple method for explicitly constructing the 
N + 1 sets of MUB. 

The association between striations and MUB also ap- 
pears in a transparent manner: As we mentioned above, 
each striation should be associated with a set of N — 1 
nontrivial commuting translations. As we have a total 
of TV 2 — 1 nontrivial translation operators (we do not 
include the identity) we can split them into N + 1 dis- 
joint sets containing N — 1 operators each. In fact, this 
was done explicitly above. All operators within each set 
commute and the operators belonging to different sets 
are orthogonal. When these conditions are met a power- 
ful theorem proved by Bandhyopadyay et al. holds |25j : 
If a set of N 2 — 1 traceless operators can be split into 
N + 1 orthogonal sets of commuting operators then the 
eigenbases associated with these sets are mutually unbi- 
ased. 

It is very important to notice that condition Q also 
implies that the mapping between field elements q,p and 
binary strings q, p cannot be arbitrary. Let us consider 
the ray formed by the points (so, sb) (which satisfy the 
equation bq + ap = 0). As described above, the points 
in this ray are associated with the translation opera- 
tors r(aM J ,bM 3 ). Consider in particular the horizon- 
tal ray p — 0, which contains all points (sa, 0). Apart 
from the origin (0, 0) all the points in this ray can be 
ordered according to the powers of the generating ele- 
ment uj since they can be expressed as (u; J , 0). There- 
fore, the binary strings corresponding to these points 
should be of the form (aM 3 ,0), and the origin. Choos- 
ing a=l = (l,0,...,0) the matrix M determine the re- 
maining binary strings associated with the position axis 
of phase space. These binary strings are nothing but 
the coefficient of the expansion of the field clement q in 
the canonical basis formed by the first n powers of lu 
(of course, choosing a different binary string a for the 
first element is tantamount to a change of basis). The 
binary strings associated with the momentum axis can 
also be determined in this way. Consider the vertical ray 
q = which contains the origin and all the points of the 
form (0, lu j ). These points are associated with the binary 
strings (0, bM 3 ). and the origin. Again, the choice of the 
binary string associated with the first point in the ray is 
arbitrary and we can choose it to be b = 1. The subse- 
quent binary strings associated with the momentum axis 
are therefore determined as p = 1M J . 

It is interesting to notice that the matrices M and M 
have dual roles: While the powers of M can be used 
to find the coordinates of field elements in the canonical 
basis, the powers of M enable us to find the coordinates of 
field elements in the dual basis. Therefore, our previous 
argument implies that the binary strings associated with 
the momentum axis of phase space are the components of 
the field elements in a basis which is a multiple of the dual 
of the canonical basis (the freedom in choosing b implies 



that the basis that should be used for the momentum 
axis is a multiple of the dual and not simply the dual). 

As we described above, the condition J7J is crucial 
in establishing a relation between a basis set in Hilbcrt 
space and a striation of phase space. However, it does not 
tell us what specific state one should associate with each 
line of the striation. In fact, this association is entirely 
arbitrary. To completely define the phase space struc- 
ture one should provide such relation between states and 
lines. In this way, using the terminology employed by 
Wootters et al. [llj one defines a "quantum net" . Thus, 
there are many quantum nets allowed by equation 0. 

One can calculate the number of allowed quantum nets 
as follows. Let us assume, for simplicity, that the asso- 
ciation between states and vertical or horizontal lines is 
fixed. There are N — 1 remaining striations and each 
of the corresponding rays can be associated with any 
of the N states of the corresponding MUB. Therefore, 
there are N N_1 quantum nets consistent with the above 
constraints and each one will lead to a different Wigner 
function as described below. 



E. Fixing the quantum net 

To specify a quantum net we need to make explicit 
the connection between every line in a striation and a 
state in a basis of the Hilbert space. For this purpose 
we only need to assign a state to the ray of the striation. 
Thus, once we do that we can obtain the states associated 
with all the other lines of the striation by applying the 
translation operator that maps the ray to the line. We 
will adopt a labeling scheme for the rays as in Fig. ^ 
The ray A is the line p — u x q and we label separately 
the horizontal (p = 0) and the vertical (q = 0) rays as 
A = h, v. 

It is useful to notice that a rank one projector onto an 
eigenstate of the translation operators associated with 
a ray can be constructed as follows: The first n points 
in each ray (besides the origin) correspond to a set of 
n translation operators which can be considered as the 
generators of the set of translations associated with the 
ray. In fact any set of n non-trivial translations belonging 
to the ray could be used as generators since the rest of 
the operators are obtained from all possible products of 
the n generators. Then the generators for ray A are the 
operators 

G[ x=j) = T(lM k , lM( fc+J ')), k = 0,...,n-l, (8) 

where, again, 1 is the binary string (10 ... 0). 

A possible, and simple, way to fix the quantum net is 
by associating each ray with the state which is an eigen- 
state of all the generators with eigenvalue +1. Recall that 
the translation operators are hermitian and unitary and 
therefore have eigenvalues ±1. Therefore this projector 
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is 

n— 1 

n (A) = F n( i+G i A) )- ( 9 ) 

fc=0 

The projectors on the other lines of the striation are ob- 
tained by translating the above one. In this way we have 
arbitrarily fixed the quantum net. Notice that the as- 
sociation of a given state to each ray is an independent 
process. Instead, once the state associated with a ray 
is chosen, the states to be assigned to the rest of the 
striation are completely determined by the covariance re- 
quirement under translations imposed in Eq. (JJJ). 

However, there operations which are not translations 
but map lines into other lines, thus leaving the phase 
space invariant: these are the unit determinant linear 
transformations. In the usual continuous case these are 
the symplectic transformations, and the usual Wigner 
function is also covariant under these more general phase 
space maps. It is then natural to try to reduce the num- 
ber of possible quantum nets by imposing additional sym- 
metries. Unfortunately, as discussed in |ll| . it is not pos- 
sible to find a faithful unitary representation for arbitrary 
linear, unit determinant transformation. However, there 
is an important transformation for which this is indeed 
possible: the rescaling (or squeezing) it w defined as 

u u {q,p) = (wg,w -1 p). (10) 

In the Appendix we describe a general quantum circuit, 
made out of simple swaps and controlled nots, defining 
the unitary operator C7 W that represents the transforma- 
tion Uuj in Hubert space. This operator is such that 

u^x a ui = x aM u^z h ul = z™' 1 . (11) 

The operator u u maps rays into rays: while leaving the 
vertical and horizontal rays unchanged it cycles through 
all the diagonal ones. In fact, using the notation intro- 
duced above, for every ray that is not vertical nor hori- 
zontal the squeezing transformation is such that u u \ — 
A — 2 (mod N — 1). The existence of this transforma- 
tion can be used to reduce the arbitrariness in defining 
the quantum net. Thus, once we fix the state to be as- 
sociated to a diagonal ray, say Pq X ~°\ then the states 
associated with all the other rays are not independent 
any more. In fact, they are fixed by the covariance re- 
quirement 

P(u u X) = U u P(X)Ul (12) 

If this condition is imposed, the arbitrariness in the quan- 
tum net is greatly reduced from N N ~ X to N. 

The operator XJ^ has another very useful property. As 
discussed above, provides the change of basis from 
the one associated with the striation A onto the one cor- 
responding to the striation A — 2 (mod N—l). Suppose 
that we are interested in explicitly obtaining the states of 
all N — 1 MUB associated with the striations which are 



not vertical nor horizontal. For this we would only need 
to find the states associated with one striation, say the 
one corresponding to A = whose generators are of the 
form = T(lM k ,lM k ). Once we find the common 
eigenstates of these generators we can easily construct 
the states of all the remaining N — 2 MUB by simply 
applying the operator U u . The simplicity of the circuit 
representing U u makes this a very efficient method to 
obtain the states of all MUB. 



F. Example 1: Phase space for n = 2 qubits 

We will show, as an example, how to construct the 
discrete phase space for systems of 2 qubits (N = 4) 
using the procedure described above. This case was al- 
ready discussed in detail by Wootters in 0, ^| but we 
include it here for completeness. The field GF(2 2 ) has 
four elements. The primitive polynomial can be chosen 
as ir(x) = x 2 + x + 1. The primitive element to is the 
root of such polynomial (which clearly does not belong 
to GF(2)). The field has then four elements: the zero 
and the three powers of u, GF(2 2 ) — {0, l,u>,l + uj}. 
The field has a canonical basis = {l,w} whose dual 
is e-i = {1 + u>, 1}. As we mentioned above, the basis fi 
to expand the momentum can be chosen as a multiple of 
£j. It is convenient to chose fi — ioei since in this case 
both bases turn out to be the same, i.e. /, = Cj. The 
phase space coordinates are the ones seen in Fig. |3 The 
generators of the horizontal, vertical and A = striation 
are: 

G { V) = T(00,10), G[ v) =T(00,01), 
G { Q h) = T(10,00), G[ h) =T(01,00), 
G { ° } = T(10,10), G { ° } =T(01,01). 

To the vertical, horizontal and A = ray we associate 
the eigenstate with eigenvalue +1 of the corresponding 
generators. The states associated with the two other rays 
A = 1, 2 can be obtained from the above by applying the 
operator U u . For the system of two qubits this operator 
is represented by the quantum circuit shown in Fig. [31 
As mentioned above, if we apply this operator to the 
eigenstates of Gq ' = Yq and G^ ' = Y\ we obtain the 
elements of the MUB corresponding to the striation A = 
2. Acting once more with this operator we obtain the 
basis corresponding to A = 1. These bases are known as 
"Belle" and "Beau" [l(j and this is a rather simple way 
to determine their states. 



G. Example 2: Phase space for n — 3 qubits 

For n — 3 position and momentum labels take values 
in the field GF(2 3 ). As mentioned above, the primi- 
tive polynomial can be taken to be ir(x) — x 3 + x 2 + 1 
and the field consists of the following elements GF(2 3 ) = 
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FIG. 2: Phase space for a system of two qubits. Two la- 
belling^ appear in the axes: One labeling associates position 
and momentum with elements of the field GF(4). The other 
labeling associates them with binary strings (which can be 
directly mapped onto quantum states). 




FIG. 3: Quantum circuit representing the operator U w for 
two qubits. This operator acts on the states associated with 
vertical lines by "moving them to the right" (i.e. mapping 
the state associated with the line q onto the state associated 
with the line u>q). Similarly, this operator acts on states as- 
sociated with horizontal lines by "moving them downwards" . 
The states associated with the lines q = and p = remain 
invariant. 



{0, 1,uj,uj 2 , 1 + uj 2 , 1 + uj + uj 2 , 1 + U),UJ + uj 2 }. The 



= {1, 



uj, us } has a dual given by 



uj 3 ei 



canonical basis ej 

gj = (uj 4 ,uj 3 ,lu 5 ). If we consider the basis fa 
we find that fa = {l,w 6 ,w}. It is interesting to notice 
that in this case it is not possible to use fa — as we did 
in the simpler n = 2 case. The phase space coordinates 
are shown in Fig. ^ The generators for the horizontal, 
vertical and A = rays are (we are using the ordering 
given in the Appendix in Table 



G, 
G 
G 
G 
G 
G 



= T(100,000), G[ h ^ — T(010, 000), 

= T(001,000), 

= T(000,100), G^ — T(000, 001), 

= T(000,011), 

= T(100, 100), Gf ] = T(010,001), 

= T(001,011). 



Again, the states associated with these rays can be ex- 
plicitly constructed using expressions like © ■ The states 
corresponding to the other rays are obtained by applying 
the operator U^, that is now represented by the quantum 
circuit shown in Fig. [S] 
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FIG. 4: Phase space for 3-qubits systems. 




FIG. 5: Quantum circuit representing the operator 
three qubits. 



for 



III. DISCRETE WIGNER FUNCTION FOR A 
SYSTEM OF n QUBITS 

We will define here the discrete Wigner function using 
the phase space structure that we introduced in the previ- 
ous section. The Wigner function provides a representa- 
tion of any quantum state in phase space, The definition 
we use will be such that the discrete Wigner function 
has the same three crucial properties than its continuous 
counterpart: PI) The Wigner function W(q,p) is real 
valued. P2) The Wigner function provides a complete 
description of the state and is such that the inner prod- 
uct between any pair of states p\ and P2 can be obtained 
as Tr(p lP2 ) = NJ2 q , p Wi(q,p)W 2 (q,p). P3) The sum 
of values of the Wigner function along any line in phase 
space is equal to the probability of detecting the state 
associated with the line. 

It is interesting to notice that once we define a phase 
space structure (or a "quantum net") the Wigner func- 
tion is uniquely determined by the condition P3). This 
can be seen as follows: Consider the point a = (q,p) and 
all the lines A that cross this point. These lines only in- 
tersect once, precisely at the point a. Therefore, the sum 
of the values of the Wigner function over all these lines 
is 

W(0)=NW(a) + J2w(P). (13) 

The condition P3) implies obviously that the sum of val- 
ues of the Wigner function over all phase space should be 
unity. On the other hand, this condition also implies that 
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the left hand side of (|13|) should be equal to the sum of 
expectation values of the projection operators associated 
with the lines A. Therefore, we obtain that P3) implies 
that the Wigner function at any phase space point is 

W («) = ^( E Tr(pP(A))-l). (14) 

A/ ctGA 

This equation can be considered to be the definition of 
the discrete Wigner function. Equivalently, the Wigner 
function can be seen to be the expectation value of a 
"phase space point operator" A(a): 

W{a) = Tr(pA(a)). (15) 

According to l|14|) the phase space point operator should 
be defined as: 

A ^ = ^( E p W-V- (16) 

A/ qGA 

The phase space point operators have several noticeable 
properties: From their definition it is clear that they are 
hermitian. Also, we can show that these operators form 
a complete basis of the space of operators. In fact, the 
Schmidt inner product between any such operators is 

Tr(A(q,p))A(q',p')) = ^S q , q ,5 p , p ,, (17) 

where 5 XtX i denotes the Kronecker-delta symbol (the va- 
lidity of this equation follows directly from 116fl using the 
properties of the phase space lines discussed above [ill). 

Expressing the Wigner function in terms of the phase 
space point operators is useful since it makes clear the 
reason why properties Pl)-P3) are satisfied: PI) is a con- 
sequence of the fact that A(q,p) are hermitian operators. 
P2) follows from the fact that A(q,p) form a complete 
orthonormal basis of the space of operators. P3) follows 
from the fact that the sum of A(q,p) over a line is noth- 
ing but the projection operator onto the state associated 
with that line. 

Using the fundamental equation J7J it is clear that the 
phase space point operators can be obtained as a trans- 
lation of the operator associated with the phase space 
origin as 

A(a) = T a A(Q)Tl (18) 

In turn, the phase space point operator A(0) is obtained 
as a sum over projection operators onto the states asso- 
ciated with the rays (lines crossing the origin) : 

A (°) = ^( E P o A) -I)- (19) 

A/ OGA 

An explicit expression of ^4(0) as a function of the trans- 
lation operators can be obtained by using @ in (|19|l . For 



n = 2 qubits this operator can be written as: 

A (°) = tq{ Epm+^d+ 

i,j=0 

2 

+ U " l l Y o + Yi+ Y a Y{\Ul m - I}. 

m=0 

For n — 3 qubits the form of A(0) is more involved and 
turns out to be given by 

A (°) E (x' xixi+ziz{zD+Y^ 

i,i,fe=0 m=0 

U™[ Y 1 + (X X Z 2 + Z X Y 2 - Y 1 X 2 )(I + Y )]Ul m - I}. 

Finally, it is worth stressing that as the operators A{a) 
form a complete basis, one can expand the density matrix 
in such basis and obtain 

p = Nj2w(a)A(a). (20) 

a 

Thus, the Wigner function W(a) is nothing but the co- 
efficient of the expansion of the state p in the basis of 
phase space point operators A(a). 

In the following sections we will display the Wigner 
function of various quantum states. However, it is clear 
that quantum states corresponding to lines have Wigner 
functions with simple properties. Thus, using (|16|) we 
can show that the Wigner function of the quantum state 
P(A) is equal to 1/N over all points in the line A and is 
equal to zero elsewhere. 

Wigner functions are useful to compute expectation 
values of operators. In particular, computing expecta- 
tion values of translation operators turns out to be par- 
ticularly simple. Thus, we can show that 

Tr(pT f3 ) = f J2w(a)(-ir^, (21) 

a 

where 

a A (3 = q Q • P/3 - q/3 • Pa 

= ^2q a iP/3i- qpiPai, (22) 

i 

is the analogous of the vector product in phase space 
(that is equal to the area enclosed by the triangle formed 
by the origin and the points a and (3). The function fp 
depends on the point f3 and on the quantum net being 
defined as fp = NTi(A(0)Tp). It is useful to point out 
how to compute fp: Let us denote the projector onto 
the state associated with the ray that crosses the point 
(3 as Pxp- This state, which is fixed by the choice of 
quantum net, is an eigenstate of Tp with eigenvalue given 
by ffj which, therefore, is always equal to ±1 and can be 
expressed as: 

fp = Tr(TpP x? ). (23) 
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There is another identity between Wigner functions 
that turns out to be useful in some calculations. For pure 
states the following identity is valid for any translation 
T a (or any other operator): 

|Tr(pT a )| 2 =Ti(pT aP Tl). (24) 

This identity, when written in terms of Wigner functions 
reads 

| £ W{p)(~lT^\ 2 = £ W(I3)W((3 + a). (25) 

This equation provides a necessary condition for the 
Wigner function to describe a pure state. Although the 
condition is not sufficient, it is useful because of its sim- 
plicity. In some applications this, together with symme- 
try arguments turns out to be enough to determine the 
value of W(a) in all phase space points. A necessary and 
sufficient condition for the Wigner function to describe a 
pure state is obtained by writing the equation p 2 = p in 
terms of Wigner functions. This implies, 

W(a) = N 2 Y, T a ^W(l3)W(j), (26) 
0i 

where V a p 1 = Tr(A a ApA 7 ). Computing these three- 
point coefficients is rather involved. Therefore, imposing 
the above necessary and sufficient condition is not a prac- 
tical way to proceed to find constraints on the possible 
values the Wigner function. 

The discrete Wigner function described in this section 
has many properties that are similar to its continuous 
counterpart. However, there are also some differences 
that are worth pointing out. For the continuous Wigner 
function (and also for other versions of discrete Wigner 
functions Q) the phase space point operators are both 
hermitian and unitary operators (up to a normalization). 
In fact, in such case the phase space point operator is a 
displaced reflection operator. For this reason, the Wigner 
function at any given phase space point can be measured 
by using an interesting tomographic technique that is 
described in some detail in [2 (and was generalized in 
|17|'). This tomographic method is a simple application 
of a more general technique to determine the expectation 
value of a unitary operator. The method does not re- 
quire to detect all marginal distributions (i.e. to perform 
a complete tomographic determination of the quantum 
state). However, the discrete Wigner function we use in 
this paper does not have this property. Thus, phase space 
point operators are generally obtained by displacing the 
operator A(0), which is not unitary. For this reason, it 
is not possible to directly measure the Wigner function 
at any phase space point using the tomographic scheme 
described in |2j. The way to determine the value of W(a) 
is by using equation l|14l) : W(a) is fixed once we know 
the probabilities for all states associated with the lines 
that contain the point a. 

There is another difference between ordinary Wigner 
functions and the ones we are describing here. The rela- 



tion between phase space point operators A(a) and trans- 
lation operators Tp is usually given in terms of a Fourier 
transform (see, for example, @|)- In our case the relation 
is by means of a different transformation, which is related 
to the Hadamard transform: 

Tp = f J2 (-l) aAf3 A(a). (27) 

a 

As mentioned above, the factor fp = NTr(TpA(0)) de- 
pends on the quantum net and takes values which are 
equal to ±1. Thus, translations and phase space point 
operators relate to each other by means of an Hadamard- 
like transform. It is also interesting to notice that the 
exponent a A j3 can be written in terms of GF(2 n ) in- 
variant objects as follows. This exponent is defined in 
terms of the binary strings defining position and momen- 
tum coordinates of the phase space points as a A (3 — 
■ P/3 — <3/3 • Pq- As we mentioned above, the basis that 
we use to order the momentum axis should be a multi- 
ple of the dual of the canonical basis. Let this basis be 
fi = us &i, i.e. the power k indicates what multiple of 
the dual basis fi is. Then, the above exponent can be 
shown to be identical to 

aAf3 = Tr(Lu- k (q a pp~qp Pa )). (28) 

The right hand side of the above equation is a basis in- 
dependent expression which is entirely written in terms 
of field operations. 

Finally, it is also worth noticing that equation (|27|l can 
be inverted and the phase space point operator can be 
written in terms of translations as 

A ( a ) = ^E(- 1 ) QA/3 ^ T ^ ( 29 ) 



All the dependence of the operators A(a) on the quantum 
net is contained in the function fp. Thus, this function 
which, as mentioned above takes values which are ±1, 
entirely defines the quantum net. 

IV. CONSTRUCTING THE WIGNER 
FUNCTION FROM THE STATE SYMMETRIES 

Having constructed the phase space representation for 
systems of qubits a natural question arises: For what 
kind of problems one expects this to be a useful tool? A 
first attempt to answer this question will be presented in 
this Section. Thus, we expect this tool to be of some use- 
fulness when the quantum state and/or the evolution op- 
erator have some degree of symmetry under phase space 
translations. In this Section we will discuss three spe- 
cific examples for which this is indeed the case. First, 
we will analyze the phase space representation of some 
maximally entangled states which are defined precisely 
as eigenstates of translation operators. As a first exam- 
ple we will analyze the case of Bell states for a system 
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of n — 2 qubits and show that these states have rather 
simple Wigner functions. A more interesting example is 
the case of stabilizer error correcting codes |2J|. These 
codes are also naturally formulated in terms of transla- 
tion operators. Indeed, the code space is defined as the 
space of common eigenstates of n — 1 commuting opera- 
tors Sj (j = 1, . . . , n — 1). These operators, which define 
the stabilizer of the code, are nothing but phase space 
translation operators. Below, we will discuss the phase 
space representation of the simplest such code (the three 
qubit error correcting code against phase errors). In this 
case, as we will see below, the usefulness of the Wigner 
representation turns out to be less obvious. Finally, we 
present the phase space representation of the celebrated 
'Mean king problem' [2l],|23]. This problem has a rather 
appealing solution when formulated in phase space. 

A. Bell states 

The Bell basis for a system of two qubits is formed by 
the states 

|*±) - -±=(|10>±|01», 

|*±) - -^(|00)±|11». (30) 

It is simple to show that these states are common eigen- 
states of the translation operators Tr u 2 q\ = XqXi and 
T(vi,ui 2 ) = ZqZ\ (the eigenvalues of both operators are 
±1). For this reason, one expect these states to have a 
simple phase space representation. Indeed, the Wigner 
function of these states must be invariant under the 
translations Ta^q) and T(q j0J 2\. The action of these op- 
erators is simple: T( u a q\ interchanges the vertical lines 
p = and p = uj 2 (and also interchanges the lines p = 1 
and p = uj). Similarly T^^ 2 ) interchanges the corre- 
sponding horizontal lines. The total number of points in 
phase space for two qubits is N 2 = 16. Each of the two 
above symmetries can be use to cut in half the number of 
independent parameters that define the Wigner function 
of a Bell state. Thus, we are left with only four parame- 
ters which characterize the Wigner function as shown in 
Fig. |SJ It is worth mentioning here that the above results 
are independent of the quantum net (the only assump- 
tion we made concerned the association of vertical and 
horizontal lines with the corresponding basis). Indeed, 
there are 4 3 quantum networks with Wigner functions 
having the symmetries shown in Fig. 

All the four Bell states have Wigner functions with 
the above symmetries. To find the Wigner function of 
each Bell state we need to impose some other conditions 
that constrain the possible values of the four parame- 
ters a, 6, c, d. There is an obvious constraint imposed by 
normalization, i.e. the sum of all values of the Wigner 
function must be unity. This implies that 

a + b + c + d=\. (31) 
4 
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FIG. 6: Wigner representation for Bell states, depends on four 
parameters. Using the symmetries of the Bell states we reduce 
our 4x4 grid of independent values to the determination of 
4 real values a,b,c and d. 

As mentioned above, Bell states are eigenstates of the 
translation operators with eigenvalues ±1. For the state 
|$ + ) both eigenvalues arc +1. Imposing this condition 
is equivalent to requiring that the expectation value of 
the translation operators is equal to +1. The expecta- 
tion values of TVw 2 ) an d ^ can be easily computed 
from the Wigner function using Q21JI (in such case, for 
both translation operators we have f@ = +1). Thus, this 
condition gives us the following two equations: 

a + b — c — d = — , 
4 

a + c-b-d =-. (32) 

Equations (|31|l and l|32l) are a simple linear system of 
three equations for four unknowns. Therefore, we still 
have a one parameter family of solutions for the Wigner 
function of each Bell state (the linear system determining 
the Wigner function for the other three Bell states |$_) 
and |\&±) is formed by (|3 1|) and two equations analogous 
to with the corresponding ± signs associated to the 
eigenvalues) . 

It is interesting to notice that the Wigner function of 
Bell states can be further constrained by imposing the 
following condition: Any Bell state is mapped onto an or- 
thogonal state by a unitary operator that anti-commutes 
with either Tjuo; 3 ) or Tr u i m. This is indeed the case for 
the translations given by Xq, X\, Zq, Z\. In fact, apply- 
ing the translation X$ maps \^±) onto |$±). Imposing 
that the Wigner function of the state translated by X$ is 
orthogonal to the one associated with the original state 
(i.e., that ^2 a W'(a)W(a) — 0) is equivalent to the fol- 
lowing equation: 

ab + cd = 0. (33) 

This equation, together with the above linear set, give 
two possible solutions: {a — j,b = c = d = 0} and {a = 
b = c = —d = I}. The two possible Wigner functions are 
shown in Fig. |7| The Wigner functions associated with 
the other three Bell states are obtained from the above 
by applying the corresponding translation operators (for 
example, the one for |Vi?+) is obtained by applying Xq, 
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FIG. 7: The two possible Wigner representations of the Bell 
state |3>+) = "^7f(|00) + |11)). These representations were 
constructed using the symmetries of the state. The other 
Bell states can be constructed by applying the translations 
Xi, Z x and XiZi. 



which corresponds to interchanging the first column with 
the second one and the third with the fourth). The two 
solutions presented here are the only ones allowed by the 
4 3 possible quantum nets. 



B. Phase space representation of a quantum error 
correction code 



Quantum error correction codes have been under in- 
tense investigation during recent years (2?|]. We will not 
present a detailed introduction to the theory of quantum 
error correction but simply introduce the necessary ele- 
ments to study their phase space representation. A wide 
class of such codes is defined in terms of a stabilizer as 
follows: Let us consider the simplest case in which we en- 
code one qubit of quantum information using n physical 
qubits. The space of encoded states is a two-dimensional 
subspace of the total Hilbert space formed by the set of 
common eigenstates of the operators Sj, j = 1, . . . , n — 1. 
Such operators are, in the jargon defined above, a set 
of commuting translation operators since they are tensor 
products of Pauli operators acting on each qubit. The 
stabilizer is chosen in such a way that the code corrects a 
set of errors Ei which are also translation operators. The 
code will correct against errors Ei, i — 1, . . . , 2™~ 1 — 1 
if the encoded states \4>l) are mapped by the errors Ei 
onto subspaces which are mutually orthogonal for differ- 
ent values of i. 

The use of phase space representation in this con- 
text seems to be natural. Indeed, as encoded states 
are eigenstates of translation operators their correspond- 
ing Wigner functions must be invariant under the same 
translations. Moreover, as errors are also translation op- 
erators, when translating the Wigner function of an en- 
coded state by a correctable error one should obtain a 
Wigner function which is orthogonal to the original one. 
Therefore, one expects (perhaps naively) the phase space 
method discussed above to provide some insight into er- 
ror correction. Also, the task of finding the Wigner func- 
tion of encoded states could be accomplished using the 



same ideas we described above to determine the Wigner 
function of Bell states (i.e., using the symmetries of the 
state in the first place). 

Here we will discuss a specific example, the simplest 
stabilizer quantum error correcting code. The code en- 
codes one qubit of quantum information using n = 3 
physical qubits and corrects against errors of the form 
Zq, Z\ and Zi (i.e., phase errors). The stabilizer of the 
code is defined by the translation operators 



Si — ^6,0) — X Xi, 
$2 = ^=,0) = AiX 2 . 



(34) 



The reason why the code can correct against all Zi errors 
is that each error maps eigenstates of Si onto eigenstates 
of the same operators with different eigenvalues. Indeed, 
if we define the code space as the set of all states with 
eigenvalues +1 for both Si and S2, the action of errors on 
encoded states turns out to be the following: Zo errors 
maps encoded states onto eigenstates of Si and S2 with 
eigenvalues —1 and +1 while for Z\ and Z2 errors the 
corresponding eigenvalues are s\ = —1, S2 = —1 and 
Si = +1, S2 — —1 respectively. 

As mentioned above, the Wigner function of encoded 
states must be symmetric under the translations Si and 
£2- As was the case for Bell states, each of these sym- 
metries cuts in half the total number of independent pa- 
rameters defining the Wigner function. So, invariance 
under Si and S2 is achieved only if the Wigner function 
is identical along the following vertical lines: W(0,p) — 
W(uj 3 ,p) = W(uj 5 ,p) = W{uj 6 , p ), W(l,p) = W(u,p) = 
W(uj 2 ,p) — W(uj 4 ,p) (the identities must hold for all val- 
ues of p). Thus, this symmetry reduces the number of pa- 
rameters defining the Wigner function of encoded states 
from 8 x 8 to 8 x 2. Below, we will display the Wigner 
function of general encoded states. But before that, it 
is simpler to start by analyzing the Wigner function of 
specific encoded states. To define the states encoding 
the two logical states |0l) and we can proceed as 
follows: These states can be chosen to be encoded states 
(which are eigenstates of eigenvalue +1 of Si and S2) 
which are also eigenvalues of a third operator that com- 
mutes with the stabilizer. This operator can be chosen in 
this case as Zl = ZqZiZ2, which is also a translation op- 
erator. Imposing invariance under Zl again cuts in half 
the number of independent parameters. Thus, as Zl in- 
terchanges horizontal lines one must identify the values of 
the Wigner function along such lines (the Zl symmetry 
implies that W(q,0) = W{q,uj 3 ), W(q,l) = W(q,u 2 ), 
W(q,uj) = W(q,uj 4 ), W{q,uj 5 ) = W{q,cu 6 )). There- 
fore, the Wigner function of the two logical states can 
be parametrized with 8 parameters and has the symme- 
tries displayed in Fig. |S] 

We can apply the same reasoning we used to find the 
Wigner function of Bell states above to find extra condi- 
tions on the above eight parameters. First, we can im- 
pose the normalization condition and the fact that logical 
states have eigenvalue +1 of the translation operators Si 
and S2- For the case of the state |0l), the eigenvalue of 
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After some algebra we find that the possible solutions 
must obey the following conditions: 



000 100 010 001 101 in no on 

FIG. 8: Wigner function for a logical state. The number of 
independent variables in this function can be reduced by using 
the symmetries of the state under the stabilizer translations. 
In the 8x8 grid there are only 8 independent variables to be 
determined (a, b, c, d, e, f and g). 



Zl, is also equal to +1 and the corresponding conditions 
turn out to be 

a + b + c + d + e + f + g + h = -, 

8 

a + b + c + d — e — f — g — h = — , 

8 

a + b— c— d + e + f — g — h = — , 

8 

a-b + c-d + e-f + g-h = ~. (35) 

8 

The conditions defining the logical state |1) ^ are the same 
as above except for the last equation where the right hand 
side is —1/8. It is worth noticing that once we have the 
Wigner function of the logical state |0l) the one for |lj,) 
is obtained by translating it with the operator Xq which 
interchanges vertical lines. 

The above linear system still allows for a four parame- 
ter family of solutions. Extra conditions can be imposed 
in the same way as we did for Bell states. Indeed, er- 
rors Z , Zi and Z 2 are such that they transform encoded 
states into orthogonal states. Imposing the orthogonality 
with the translated Wigner functions yields the following 
set of equations 

ae + bf + eg + dh = 0, 
ac + bd + eg + fh = 0, 
ag + bh + ce + df = 0. (36) 

Finally, to find solutions that correspond to quantum 
states we should impose the condition that the sum of 
values of Wigner function along arbitrary lines should 
always be non-negative. Also, we impose the condition 
for the state to be pure (i.e., N J] a W 2 (a) = 1), which 
is equivalent to 



b = — — a, f = — e, d = — c, h = —g. 



(38) 



We find eight solutions with all the desired properties. 
They are 

a = — , e = c = g = 0; 



c=^e=g=Q 



q = — , c — 
y 16' 


e = 
1 




c = e = g = 


32' 




3 




1 


— ,c — e — 
32' 


-g = 


32' 


3 




1 


— , c = — e 
32' 


= 9 = 


32' 


3 




1 


32'~ C = e 


= 9 = 


32' 



(39) 



So far we did not impose any condition on the quantum 
net. In fact, some of the Wigner functions obtained in 
this way do not have the property of covariance under the 
operation U u . Imposing this condition we are left with 
only four solutions which are the ones corresponding to 
the last four equations. In Fig. El we show one of these 
solutions, and the solution where the Wigner function 
only takes positive values. 

Finally, we completely specified the quantum net by 
choosing the state corresponding to the main diagonal 
to be Zi|Ao), where |Ao) is the eigenstate with eigenvalue 
+1 of the three generators corresponding to this line. For 
this quantum net the Wigner function of the encoded 
state |0)i is the one shown at the top of Fig. [5] For 
this case, we also obtained the Wigner function of the 
most general encoded state \<f>L) = oi\0l) + This 
Wigner function is displayed in Fig. ^]and is completely 
determined by the following four functions: 



f2(a,0) 



a 



32 1 

>' 2 



3|/3| 2 + (2 + i)a/3* + (2 - i)a* 0\, 
\(3\ 2 + i(ap* -a*0)], 
\p\ 2 -i(af3*-a*P)}, 
■ 3|/3| 2 - (2 + i)ap* - (2 - i)a*0\. 



c 2 +d 2 + e 2 + f + g 2 + h 2 



1 

64' 



(37) 



The Wigner function shown in Fig. EHhas all the sym- 
metries defining encoded states (i.e., it is invariant under 
interchange of vertical lines corresponding to the transla- 
tion operators Si and S% and is mapped onto an orthogo- 
nal state when translated by errors Z , Z\ and Z%). Our 
naive expectation was that these properties were going 
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FIG. 9: Two possible Wigner representations of the state |0)l . 
The one at the top, where {a = c = e = g= -L}, is covariant 
under the action of t7 w . The quantum net is defined by asso- 
ciating the eigenstate Zi|Ao) to the main diagonal line (where 
| Ao) is the eigenstate with eigenvalue +1 of the generators). 
The one at the bottom, where {a — g, e = c = g — 0}, is not 
covariant under U u . The color convention is such that black 
(white) regions correspond to positive (negative) values of the 
Wigner function. 

to be more evident in the solution we obtained. How- 
ever, this is not the case, which casts doubts about the 
usefulness of the phase space representation for quantum 
error correction. 



C. The Mean King Problem: a phase space 
solution 

Here we present a third application of phase space tools 
to quantum information. We will solve the so-called 
mean king problem, that was first presented in |2lJ ■ The 
formulation of the problem is the following: A physicist 
must prepare a spin 1/2 particle in a state of his choice. 
Then he should give the particle to the mean king. The 
king makes a projective measurement of one of the three 
Cartesian components of the spin of the particle (i.e., the 
king measures either the X, the Y or the Z observable). 
Then, the king gives the particle back to the physicist 
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FIG. 10: The value of the Wigner function of a general en- 
coded state a|0) l +/3|1) l in the first two columns of the phase 
space. The quantum net is the same one used at the top of 
Fig. E] The other columns can be obtained from the invari- 
ance of the state under the stabilizer translations Si and S2 ■ 



who is allowed to perform any operation on it. Finally, 
the king announces what observable was measured in his 
laboratory. The physicist would only save his life if he 
is able to retrodict the result of the king measurement 
based on the results of the measurements performed in 
his own laboratory before knowing the observable mea- 
sured by the king. This problem was also extended to 
cases where the physicist is given quantum systems with 
a Hilbert space with dimension which is a power of a 
prime (22]. The solution of the mean king problem is 
only possible if the physicist entangles the particle to be 
sent to the king with another identical particle he keeps 
in his own laboratory. 

We will analyze this problem using phase space meth- 
ods, which seem to be well suited for this purpose. In- 
deed, maximally entangled (Bell) states are naturally 
represented in phase space as discussed above. Moreover, 
the observables measured by the king are also naturally 
represented in phase space since they are translation op- 
erators. For this reason, one may suspect that phase 
space methods may enable a simple solution to the prob- 
lem. This is indeed the case, as we will now discuss. We 
will not present here the usual solution to the problem 
but will attempt to present a solution entirely based on 
phase space. Let us consider the initial state prepared by 
the physicist to be the Bell state |$+) whose Wigner func- 
tion was displayed above. When the king measures one 
of the three components of the spin of the first particle he 
can obtain one of two values (±1/2). Each of these mea- 
surements can be viewed as the projection of the original 
state onto a basis of "line states" : Vertical lines are as- 
sociated with the measurement of Z, horizontal lines are 
associated with the measurement of X and the striation 
containing the ray corresponding to the main diagonal of 
phase space corresponds to the measurement of Y. From 
the Wigner function of the above Bell state it is clear 
that the initial state has non-vanishing projection only 
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on two states of each of these three striations. Therefore, 
as a result of his measurement, the king prepares one of 
six states which are displayed in Fig. El The quantum 
net we use here is such that the states corresponding to 
the two vertical lines are (right vertical line denoted 
as v%) and |00) z (left vertical line, denoted from now on 
as the two states corresponding to the horizontal 
lines are \ll) x (top horizontal line: h%) and \00) x (bot- 
tom horizontal line, denoted as h 2 ) and the two states 
corresponding to the diagonal lines are |10) y (diagonal 
not crossing the origin: di) and |01) y (diagonal crossing 
the origin, denoted as d 2 ). 




11 




11 



FIG. 11: A Wigner representation of the state |<£+). At the 
right we show the possible states that the king could give us 
according to his measurement. The six states are associated 
to the vertical lines vi and V2, to the horizontal lines hi and 
fi2 or to the diagonal lines di and d,2. 



Once the king makes his measurement and gives the 
particle back to the physicist, he should device a mea- 
surement scheme which would enable him to retrodict 
the result of the king's measurement once he is informed 
about the measured observable. What can the physicist 
do? The solution is to measure a collective observable 
of the two particle system with four distinct eigenvalues. 
Each of the four eigenstates of this observable must have 
vanishing overlap with one and only one of the two states 
generated by each measurement. In such case, when the 
king announces what was the measured observable the 
physicist will always be able to infer what was the mea- 
sured result (since there is only one result of every king 
measurement which would be consistent with the result 
of his own measurement). For example, one such state, 
which we denote as should be orthogonal to the line 
states associated with the upper horizontal line (hi), to 
the rightmost vertical line (vi) and to the diagonal line 
not crossing the origin (d\). As mentioned above these 
three line states are |ll) z , |H)x and 1 10)^. Below, we will 
show how to obtain the Wigner function of this state. 
It is important to notice that once we obtain this func- 
tion we can find the Wigner functions associated to the 
other three states that complete the basis of the Hilbert 
space by implementing simple translations. This is the 



case for the following reason: The six line states rele- 
vant for our construction are connected by phase space 
translations. For example, the operator T^ qj = XqXi 
interchanges the two vertical states and the two diago- 
nal states. Similarly, the operator T^q.^^ = Z^Z\ inter- 
changes the two horizontal states and the two diagonal 
states while ^2^2) interchanges the two vertical states 
and the two horizontal ones. Therefore, it is natural to 
require that |cpi) and its three translated descendants 
form an orthogonal basis of the Hilbert space. This basis 
defines the observable to be measured by the physicist. 

To find the Wigner function of the state \<px) we need 
to determine its value in all the 4x4 points of the phase 
space grid. For this purpose, we can proceed much in 
the same way we did in the previous sub-sections. Thus, 
we impose the following conditions: The first condition 
tells us that the sum of the values of the Wigner function 
along the three lines hi, v\ and d\ must be equal to zero. 
Fig. IT^l shows the three lines that define this states. This 
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FIG. 12: The Wigner function of the state \ipi) that defines 
the solution to the mean king problem should add up to zero 
along the lines indicated in the Figure. As it is evident from 
the graph, is natural to impose that the Wigner function is 
invariant under reflection about the main diagonal. 

condition imply the following set of three linear equations 



W(a) = w i a ) = E w i a ) = °- ( 40 ) 



The second condition arises by imposing that the three 
translation operators T( w 2 ), T(q,u) 2 ) and T^^j map the 
state into an orthogonal one. This is equivalent to impose 
that the expectation value of these operators is equal to 
zero. In turn, this condition is naturally expressed in 
terms of the Wigner function as: 



(X 1 X 2 ) = 
(ZiZ 2 ) =0 

(X1Y2) = 



E 



W(a) 

E w ^) - 



E ^w. 

a^hi,h 2 

W(a), 



E w (< 



E 
E 



W(a) 



It is simple to show that these conditions, together with 
the normalization condition, imply that the sum of the 



Wigner function along the lines h-2, Vi and di must be 
equal to 1/2, i.e. 

J2 W(a) = w ( a ) = E = \- ( 41 ) 

a£/i2 a£v2 a£d,2 

Also, for the same reason the sum of the Wigner function 
along the remaining two lines of each striation (i.e., ^3,4, 
and g?3,4) must be equal to 1/2, i.e. 

J2 W(a)= ]T W(a)= W ( a ) = \ 

a£h,3,h.4 a£v3,V4 (x^d^.d^ 

(42) 

Clearly, the above conditions are not enough to com- 
pletely determine the Wigner function. However, we can 
obtain a solution by noticing that the state \tp\) has an 
extra symmetry. In fact, it is invariant under the opera- 
tion P12H1H2 which interchanges the two particles after 
applying a Hadamard transformation to each of them. 
This is so because the Hadamard operation applied to 
both qubits interchanges lines hi with v\ and d\ with 
d.2- On the other hand, the permutation leaves hi and v\ 
invariant while interchanging d\ and £?2 . Using this sym- 
metry, together with the previous results, we can show 
that the sum of the Wigner function along lines /13 and /14 
has the same value (similarly for the remaining vertical 
and diagonal lines). 

So, all the conditions discussed so far can be summa- 
rized as follows: In each of the three striations (vertical, 
horizontal and "diagonal" ) the Wigner function adds up 
to 1 /2 in one line, it adds up to 1 /4 in two other lines and 
it adds up to in the remaining line. It is simpler to find 
a solution that has the same symmetry shown in Fig. 1121 
i.e. reflection about the main diagonal of phase space. 
For this reason, the Wigner function we will obtain can 
be parametrized by ten real numbers as shown in Fig. 
ITS1 The ten parameters should obey all the conditions 
listed above, which end up giving rise to seven indepen- 
dent linear equations. To fix the solution we need to 
impose other constraints. The need for these constraints 
comes from the fact that so far we did not impose the 
Wigner function to describe a state. This can be done 
by imposing, for example, that the purity is equal to one 
(i.e., NJ2 a W 2 {a) = 1). Also, we can use the fact that 
the orthogonality between \<pi) and the translated state 
\<P2) = T(p tU 3)\(pi) implies that J2 a Wi(a)W 2 {a) = 0. 
It is worth stressing that this condition is equivalent to 
Tr(pT( iU 2)) = only if the state is pure. Finally, we can 
impose that the Wigner function can be used to compute 
the expectation value of, for example, the operator Xq in 
two equivalent ways, i.e. 

I W((3)(-l) a ^\ 2 =NJ2W((3)W(P + !)■ (43) 

where 7 = (1,0). The above set of conditions give four 
possible consistent Wigner functions. The one corre- 
sponding to the quantum net which is covariant under 
the operator U u is given in Fig. [21 
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FIG. 13: Ansatz for the Wigner function of the state \ipi) 
defining the solution of the mean king problem. It depends 
upon 10 real parameters and is symmetric upon reflection 
about the main diagonal. 
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FIG. 14: The Wigner function of the state \tpi), which defines 
the solution to the Mean King Problem. The physicist needs 
to measure an observable which is diagonal in the basis formed 
by this state and the ones obtained by translating it using 
the operators ZqZi, YqYi and XqX%. The color convention 
is such that positive (negative) values of the Wigner function 
correspond to black (white) regions. The solution is such that 
a — e — d — h = 3/16, b = c — f — 1/6, j — i = g = —1/16. 

V. CONCLUSIONS 

In this paper we reviewed the discrete Wigner func- 
tion for systems of n qubits. We showed a rather simple 
way to construct the phase space for these systems. The 
matrix representation of the finite field GF(2 n ) and the 
quantum circuit representation of unitary operators are 
two very useful tools in this context. The position and 
momentum axis of phase space are labeled with elements 
of the field GF(2 n ). To associate each vertical line with 
a computational state one associates each field element 
with the n-tuple formed by the coordinates of the field 
element in a given basis. The same criterion is applied for 
the horizontal lines. In this approach there is a unitary 
operator U u that plays a crucial role: U w is the unitary 
operator that permutes computational states in such a 
way that while fixing the state associated with the first 
vertical line (q = 0) it maps every other computational 
state onto the one located immediately "to the right" . 
The same operator maps momentum states moving them 
"downwards" and therefore corresponds to the classical 
squeezing operation that maps a point (q,p) onto an- 
other point (ujq,p/uj). It turns out that not only this 
operator can be easily represented in terms of a simple 
circuit entirely made out of control-not and swap gates. 
Also, the operator provides the change of basis between 
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the N — 1 mutually unbiased bases which are associated 
with all the (non-vertical or non-horizontal) striations. 
In fact, U u can be used to explicitly generate all the 
states of N — 2 bases given the states of a single one (for 
example, defining the states associated with the striation 
corresponding to the main diagonal of phase space one 
would obtain the remaining ones applying powers of U u ). 
It is important to stress that imposing covariance of the 
Wigner function under this kind of "squeezing" transfor- 
mation reduces substantially the ambiguity in defining 
the quantum net (but does not completely fixes it). 

We also presented here some new applications of the 
phase space formalism. We studied the phase space rep- 
resentation of quantum stabilizer codes (or states). As 
these states are eigenstates of a set of translation opera- 
tors their Wigner function exhibit some obvious symme- 
tries. The approach described in the paper, that enabled 
to find the Wigner function from the state symmetries, 
can be generalized in some obvious ways. Thus, we can 
show that for a stabilizer state the symmetry of the state 
under phase space translations implies that the Wigner 
function in the entire N x N grid can be constructed from 
the value it takes in TV points. Moreover, the value of the 
Wigner function for stabilizer states is always an integer 
multiple of 1/N 2 . Thus, for stabilizer states we can ob- 
tain a simple formula for the Wigner function as follows: 
A stabilizer state is defined to be a pure eigenstate of a 
set of N commuting translation operators. Consider the 
state ps to be such that Ti(psTp) — gp = ±1 for j3 £ S. 
That is to say, the state is an eigenstate of Tp with eigen- 
value given by gp = ±1 if the index belongs to a set S. 
For all other translation operators the same state has 
vanishing expectation value, i.e. Tr(psTp) = if j3 £ S. 
Therefore, using equation l|29[l the Wigner function of a 
stabilizer state is 



w/ s(«) = ^E^/3(-ir A/3 . 



(44) 



pes 



This result is valid for all stabilizer states and for all 
phase space points. From this equation it is evident that 
the Wigner function has the same value for all the points 
of the stabilizer: Thus, if we consider a £ S we must 
use a A (3 = in (|44|l since the corresponding transla- 
tions commute. In such case we have Ws(a £ S) = 
J2f3es ff39p/N 2 - Moreover, we can show that for every 
stabilizer S the quantum net can be chosen in such a 
way that fp — gp (for (3 £ S). Then, the Wigner func- 
tion is equal to 1/N in the points of the stabilizer and 
is equal to zero everywhere else. These features are ev- 
ident in the examples we discussed above where we an- 
alyzed Bell states and encoded states of the three qubit 
error correcting code. However, the results we obtained 
are not conclusive about the potential usefulness of the 
phase space approach as a natural tool for this problem 
(even from a purely pictorial point of view). More work 
along this line is in progress. The phase space solution 
of the "mean king problem" is also naturally formulated 
in phase space due to the central role played by mutually 



unbiased measurements in this context. 

As the Knill-Gottesman theorem states 20] , the stabi- 
lizer sates and their evolution under unitary operations 
of the Clifford group can be efficiently simulated on a 
classical computer. Therefore, one could ask if this class 
of quantum computation also induces a "classical" evo- 
lution in phase space, i.e. if it can be represented by 
a classical flow mapping phase space points into phase 
space points. For this to be true, the group of phase 
space point operators A(a) should be mapped onto itself 
under unitary transformations of the Clifford group. Un- 
fortunately, this is not the case. To see this, it is enough 
to present a counter-example and this is provided by 
the set of Bell states we studied above. Indeed, Bell 
states are stabilizer states generated by applying opera- 
tors of the Clifford group to computational states (i.e., 
a Bell state is obtained by applying a Hadamard and a 
controlled-not operator to a computational state). As 
we showed above, while computational states have posi- 
tive Wigner functions with support on the vertical lines, 
the Wigner function of Bell states have negative values 
for some quantum nets. One could ask if it is possible 
to choose a quantum net such that all stabilizer states 
have positive Wigner function. The answer to this ques- 
tion is, again, negative as the following reasoning shows: 
The quantum net is fully characterized by the function 
fp = Tr(Tp Fx*). On the other hand, a stabilizer state 
represented by a projector Ps is characterized by the 
function gp = Tr(Ps Tp). For the Wigner function to be 
positive gp must be equal to fp for all stabilizer states. 
However, it is simple to show that if these two functions 
are identical for a given stabilizer state, one can always 
construct another stabilizer state for which gp ^ fp by 
using elements of the Clifford group. Therefore, for every 
quantum net one can find some stabilizer states with neg- 
ative Wigner functions. We stress that this does not con- 
tradict the conjecture formulated by Galvao in whose 
validity would imply that states with negative Wigner 
function are necessary for exponential quantum compu- 
tational speed-up. 



APPENDIX A: REPRESENTATION OF THE 
FIELD GF(2 n ) BY n x n MATRICES 

The standard definition of the companion matrix to 
the binary polynomial 

tt(x) = r + nx + r 2 x 2 + ... + r^a;" -1 + x n 



(Al) 



/ 1 
10 



M 



\ 





r 71 r 2 • • • 



\ 




(A2) 



It has the basic property that det(M — x) = ir(x) so that 
7r(M) = 0. Thus M is a root of the polynomial, and if 
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ir(x) is primitive it can be shown [23 that M 2 1 = 1 
and that M j f 1 for all j < 2 n - 1. Thus M can be put 
in one to one correspondence with the abstract elements 
bj' J of the field. Therefore the set of powers of M gives 
a representation of the field in terms of binary matrices 
with addition and multiplication in the field given by 
ordinary (binary) matrix operations. 

The first n — 1 powers constitute a canonical basis for 
the field, and arbitrary elements in this basis can be ex- 
panded as 



.4: 



n-l 

E 

i=0 



a,- AT 



(A3) 



and are then fully represented by a binary string a 
(do, ■ ■ ■ , »n-l)- 

The action of M on field elements is then 



A' = AM 



n-i 

E 

i=0 



CLiM 



'■+1 



(A4) 



Using the fact that ir(M) = (equivalent to considering 
the multiplication of polynomials modulo n(x)) we can 
rewrite this action in matrix form on the binary string 



a' = aM. 



(A5) 



The successive powers of M acting on a reference string 
a ^ generates a cyclic ordering of the binaries with 
period 2" — 1. The opaque looking definition of the trace 
of a field element in Eq. @ becomes here the ordinary 
(binary) matrix trace. 

We are also interested in the action of M on elements 
expanded in the dual basis. This is simply given by the 
transpose M. It should be noted that this matrix also 
satisfies ir(M) = and shares all the primitive properties 
of M. When acting on a reference binary, it still cycles 
through all the non-zero binaries, but in a different order. 
Table I shows the two different orderings obtained for n — 
2, 3 and 4 qubits for the indicated primitive polynomials. 

Extensive tables of primitive binary polynomials are 
available [2!j, and for moderate values of n the above 
algorithm to generate the sequences of binaries is easily 
implemented. 

We are now interested in obtaining a unitary rep- 
resentation of this action. Consider the translations 
T(a,0) = X a defined in Eq. ©. An operator U u that 
achieves this can be constructed using only two basic 
gates: CNOT tJ and SWAPy. The operator CNOT^ 
acts on qubits i and j transforming the state \xi, Xj) into 
\xi, Xi + Xj). The operator SWAPij interchanges qubits 
i and j. The operator U u can then be shown to be 



= Y[ CNOTll J| SWAP] 



(A6) 



00 


00 


10 


10 


01 


01 


11 


11 
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n(x) 


= x 3 + X 


000 
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100 
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001 
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111 
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110 


101 


011 


010 



+ 1 




TABLE I: The binary orderings obtained in the canonical ba- 
sis and its dual for the indicated primitive polynomials. In 
both columns the reference binary has been chosen (arbitrar- 
ily) as unity. 




FIG. 15: Circuit that implements the action of a primitive 
element that is a root of the primitive polynomial ir(x) in 
Eq. <AH . The CNOT gates only act if the corresponding 
coefficient is non-zero. 



represent it as a circuit acting in the basis of eigenstates 
of Z , the computational basis. 

With the help of some circuit algebra, i. e. by commut- 
ing Pauli matrices with CNOT gates we obtain 



(A7) 



which is the same as the classical action Eq. IjASjl . To 
compute the action of U u on Z h we apply a tensor prod- 
uct of Hadamard transforms 



where 



(A8) 
(A9) 



The operator is completely determined by the coef- 
ficients of 7r(x) and its action is best understood if we 



The action of this operator is again best understood as 
a circuit. The Hadamard gates just reverse the controls 
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FIG. 16: Hadamard transform of the circuit in Fig. 1151 

in the CNOT's and the resulting circuit (in the compu- 
tational basis) is shown in Fig. ^] 

With another bit of circuit algebra we find that 
K,X b yj = X h£rl and therefore 

U u Z h Ul = Z hfy ~\ (A10) 
The action of U u on the translations is then 

U u T(q,p)Ul = ±T(qM,pM~ 1 ). (All) 



The ambiguity in the sign derives from the n/2 phases 
and can be calculated but it needs not concern us here. 

In Eq. © we defined the ray passing through 
the point (a, b) as the family of commuting operators 
T(aM J , bM J ). Acting with on this family we obtain 

U u T(aM j ,bM j )Ul = ±T(aM j+1 ,bM j - 1 ). (A12) 

This is now another family corresponding to the ray 
through (aM,bM _1 ). Thus the powers of U u cycle 
through the N — 1 "diagonal" rays, however leaving in- 
variant the horizontal and vertical ones. 
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